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^ ABSTRACT 

I Aims. Results obtained in very-high-energy (VHE; E > 100 GeV) y-ray observations performed with the H.E.S.S. telescope array are used to 
I investigate particle acceleration processes in the vicinity of the young massive stellar cluster Westerlund 1 (Wd 1). 
. . Methods. Imaging of Cherenkov light from -y-ray induced particle cascades in the Earth's atmosphere is used to search for VHE y rays irom 
^ the region around Wd 1. Possible catalogued counterparts are searched for and discussed in terms of morphology and energetics of the H.E.S.S. 
' ~~' source. 

Results. The detection of the degree-scale extended VHE y-ray source HESS J1646-458 is reported based on 45 hours of H.E.S.S. obser- 
vations performed between 2004 and 2008. The VHE y-ray source is centred on the nominal position of Wd 1 and detected with a total 
statistical significance of ~ 20 cr. The emission region clearly extends beyond the H.E.S.S. point-spread function (PSF). The differential en- 
ergy spectrum follows a power law in energy with an index of F = 2.19 ± O.OSstat ± 0.20sys and a flux normalisation at ITeV of Oo = 
(9.0+ 1.4s,a, ± l.Ssys) X 10"" TeV-'cm-^s"'. The integral flux above 0.2 TeV amounts to (5.2 ±0.9) x lO-^'cm-^s"'. 
Conclusions. Four objects coincident with HESS J1646-458 are discussed in the search of a counterpart, namely the magnetar CXOU J164710.2- 
455216, the X-ray binary 4U 1642-45, the pulsar PSR J1648-4611 and the massive stellar cluster Wd 1. In a single-source scenario, Wd 1 is 
favoured as site of VHE particle acceleration. Here, a hadronic parent population would be accelerated within the stellar cluster Beside this, there 
is evidence for a multi-source origin, where a scenario involving PSR 11648^611 could be viable to explain parts of the VHE y-ray emission of 
HESS J1646^58. 

Key words. Galaxy: open clusters and associations - Galaxy: individual objects : Westerlund 1 - gamma rays: observations - galaxies: star 
clusters 
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1. Introduction 

The long-standing question on the origin and acceleration mech- 
anisms of hadronic and leptonic Galactic cosmic rays (GCRs) 
is still not settled, despite considerable progress. The detection 
of very-high-energy (VHE) y-ray emission from shell-type su- 
pernova remnants (SNRs), e.g. Cassiopeia A, RX J 17 13-3946, 



RX J0852.0-4622, RCW 86, SN 1006 (summarised in |Hinton 



& Hofmanni 2009)1, and recently HESS J 173 1-347 (HESS 



Collaboration et al.|201 l|l and Tycho's SNR ( Acciari et al. 



lOTTj l 



supports the widely accepted idea of SNRs being acceleration 
sites of GCRs. It has been noted for many years that the Galactic 
SNR population provides sufficient energy input to sustain the 
CR flux measured at Earth. The underlying theory assumes that 
electrons and protons are injected into SNR shock fronts where 
they are accelerated via the diff' usive shock acceleration pro- 



cess up to energie s of ~10'^eV (|Krymskii||l977 



Axford et al. 



T977| |BelI||T978| [Blandford & Ostriker||1978f rThe abiUty of 
SNRs to accelerate electrons up to the so-called knee in the 
differential energy spectrum of the GCRs and our common be- 
lief that this holds for protons, too, constitute the paradigm that 
SNRs are the long-thought sources of GCRs. In interactions 
with the ambient medium, i.e. matter and electromagnetic fields, 
these GCRs then produce VHE y rays which can be detected 
by current imaging atmospheric Cherenkov telescope (lACT) 
systems, e.g. H.E.S.S., MAGIC, VERITAS or CANGAROO-III. 
Additionally, evolving SNRs could explain the chemical compo- 
sition up to the knee region. Furthermore, core-collapse super- 
novae cou ld explain observed overabund ances of some isotopes, 
e.g. ^^Ne ( [Higdon & Lingenfelter|2003 1. However, recent stud- 
ies applied to RXJ17 13-3946 highlight potential problems for 
a dominant hadronic interpretation for this object ( jEUison et al. 
12010) and motivate the search for other acceleration sites and 



processes. 

SNR shells are not the only sites in the Galaxy where GCRs 
can be produced via diffusive shock acceleration. One alterna- 
tive scenario is particle acceleration in strong shocks in colliding 
wind binaries (CWBs). Massive stars are to a large extent bound 
in binary systems (e.g. |Zinnecke"r||2003 [ |Gies||200"8| , generally 
exhibit high mass-loss rates (lO^^Mgyr"' - lO^^Mgyr"') and 
drive strong supersonic winds with velocities of the order of a 
few 10^ km s"' . When these winds collide in a stellar binary sys- 
tem they form a wind-wind interaction zone where charged par- 
ticles can be accelerated to high energies (e.g. Eichler & Usov 
|1993 i Electrons can then up-scatter stellar photons present in 
the wind collision zones via the inverse Compton (IC) process to 
GeV energies (M ucke &Poh l 2002; Manolakou et al. 2007). On 
the other hand, relativistic nucleons can inelastically scatter with 
particles in the dense wind and produce ;r°s which subsequently 
decay into VHE y rays (Benaglia & Romero H2003 Bednarek 
2005 Domingo-Santamaria & Torresl2006 ; Reimer et al.|2006 



Apart from acceleration in binaries, GCRs can be accelerated in 
the winds of single massive stars (e.g. Montmerle|1979| l. 

Another scenario involves collective stellar winds: It is com- 
monly accepted that the bulk (if not all) of the core-collapse SN 
progenitor stars and CWBs evolve from collapsing gas conden- 
sations in giant molecular clouds (e.g. |Zinnecker & Yorke 2007) l 
and mostly remain close to their birthplaces in groups of loosely 
bound associations or dense stellar clusters. When the winds of 
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multiple massive stars in such systems collide they form a col- 
lective cluster wind which drives a giant bubble (<9( 100 pc)), also 
referred to as superbubble (SB), filled with a hot (r ^ 10'' K) 
and tenuous (n < 1 cm"^) plasma (e.g. [Weaver et al.[l977[(Sihch| 



et al. 2005 1. At the wind interaction zones, e.g. at the termi- 
nation shock of the stellar cluster wind, turbulences in form of 
magneto-hydrodynamic (MHD) fluctuations and weak reflected 
shocks can build up. Unlike SNR shock fronts and CWBs where 
GCRs are accelerated through the P' order Fermi acceleration, 
turbulences in SB interiors can accelerate particles to ve ry high 
energies also via the 2""* order Fermi mechanism (e.g. [Bykov 



2001 ). Moreover, after a few million years, supernova explosions 
of massive stars (M > 8 Mq) in the thin and hot SB environment 
eventually lead to efficient particle acceleration at the boundary 
of the SB or at MHD turbulences and further amplify existing 
MHD turbulences (e.g. Ferrand & Marcowith 2010 and refer- 
ences therein). The interaction of these GCRs with the ambient 
medium including molecular clouds or electromagnetic fields 
leads to the production of VHE y rays which can then be stud- 
ied on Earth. Therefore, stellar clusters are promising targets to 
study acceleration and propagation processes of GCRs. 

One of the most prominent objects among stellar clusters 
in the Galaxy is Westerlund 1 (Wd 1). After its discovery in 
1961 (Westerlund 1961[ l subsequent observations have estab- 
lished Wd 1 as the most massive stellar clust er in our Galaxy 
placing a lower limit on its mass of 10^ M© (Clark et al.|2010| . 
An unprecedented accumulation of evolved massive stars is 
found without indication of the presence of an early-type main- 
sequence star Amongst the most massive stars, 24 Wolf-Rayet 
stars (binary fraction > 62%) have been detected and a number 
of ~150 OB super- and hypergiants (binary fraction ~ 30%) is 
expected ( Crowther et al.|2006 Dougherty et al.|2010 and ref- 
erences therein). 

The analysis of Chandra data revealed an arc minute-scale 
extended diffuse X-ray emission (Muno et al. 2006b) which is 
only seen for a few young stellar associations in the Galaxy, for 
example RCW 38 (|Wolk et al.|[20 02) and possibly the Arches 
cluster (Law & Yusef-Zadeh 2 004| l as well as in the Large 
Magellanic Cloud in 3 Doradus C ( [Bamba et"aL]|2004 ) and 
DEM LI 92 (Cooper et al. 2004). The total X-ray luminosity of 
the observed diffuse emission within 5' of Wd 1 is dominated by 
its non-thermal component and amounts to Lx ~ 3 x lO^'* erg s"' 
which represents just a fraction of 10"^ of the total mechani- 



cal power in this system (Muno et al. 2006b) . However, mod- 
els as in Oskinova ( 2005 i predict a thermal X-ray luminosity of 
~ 10^^ erg s ' for stellar clusters comparable to Wd 1, which was 
clearly not observed by Chandra for Wd 1 . As for previous ob- 
servations, there remains the open question into which channel 
most of the unobserved energy is dissipated. 

The detection of VHE y-ray emission from HESS J 1646- 
458 was initially reported in |Ohm et al. ( 2010a| l and Ohm et al 



(2010b I. This paper focuses on a detailed spectral and morpho- 
logical study of the emission region and investigates a possible 
multi-source origin. An in-depth search for plausible counter- 
parts is conducted and possible acceleration-mechanism scenar- 
ios are elaborated. 



2. H.E.S.S. Observations and Data Analysis 

Given the large ~ 5° field of view (FoV) combined with the good 
off-axis sensitivity, observations with H.E.S.S. are perfectly 
suited to cover the vicinity of Wd 1 and allow for the detailed 
morphological study of extended sources such as HESS J1646- 
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458. Thereby any large-scale non-thermal VHE y-ray emission 
around Wd 1 can be probed. 

2.1. The H.E.S.S. Experiment 

The High Energy Stereoscopic System (H.E.S.S.) is an array of 
four imaging atmospheric Cherenkov telescopes located in the 
Khomas Highland of Namibia, 1 800 m above sea level. The tele- 
scopes are identical in construction and each one is comprised 
of a 107m^ optical reflector composed of segmented spheri- 
cal mirrors. These focus the incident light into a fine-grained 
camera built of 960 photomultiplier tubes. By means of the 
imaging atmospheric Cherenkov technique (see e.g.jHillas 1985| l 
Cherenkov light, emitted by the highly-relativistic charged par- 
ticles in extensive air showers, was imaged by the mirrors onto 
the Cherenkov camera. A single shower was recorded by multi- 
ple telescopes under different viewing angles. This allowed for 
the stereoscopic reconstruction of the primaries' direction and 
energy with an average energy resolution of 15% and an event- 
by-event angular resolution better than 0.1° (Gaussian standard 
deviation, |Aharonian et al.|2006a| l . 



2.2. The Data Set 

The region around Wd 1 was observed during the H.E.S.S. 
Galactic Plane Survey (GP S) in 2004 and 2007 ( jAharonian etaL] 
2006c Chaves et al. 2008 1. Additionally, follow-up observations 
pointing in the direction of Wd 1 have been performed from May 
to August 2008. Data taken under unstable weather conditions 
or with malfunctioning hardware were excluded in the standard 
data quality selection procedure ( [Aharonian et al.|2006a) ). Also, 
pointed observations on Wd 1 at very large zenith angles of more 
than 55° were excluded due to systematic effects in the descrip- 
tion of the camera acceptance at such low altitudes for an ex- 
tended source like HESS J 1646^5 8. After quality selection and 
dead time correction the total observation time of 45. 1 hours was 
reduced to a live time of 33.8 hours. Observations have been car- 
ried out at zenith angles from 21° to 45° with a mean value of 
26° and an average pointing offset from the Wd 1 position of 
1.1°. 



rameter results in an improvement in terms of sensitivity com- 
pared to the H.E.S.S. Standard Analysis of ~ 20% and ~ 10% 



2.3. Analysis Technique 

The data set presented here was processed using the H.E.S.S. 
Standard Analysis for shower reconstruction (Aharonian et al.| 



|2006 a) and the Boosted Decision Trees (BDT) method to sup- 
press the hadronic background component (Ohm et al.„2009 ^ 
By parametrising the centre of gravity and second moments of 
the recorded extensive air shower image (! Hillas|198 5^ in multi- 
ple telescopes the shower geometry of the incident primary par- 
ticle was reconstructed stereoscopically. The directional infor- 
mation together with the measured image intensity was used to 
reconstruct the energy of the event. Since observations have been 
conducted over four years, the optical reflectivity of the H.E.S.S. 
mirrors varied and the gains of the photomultipliers changed. 
This effect has been taken into account in the spectral reconstruc- 
tion by caUbrating the energy of each event with single muon 
rings ( Aharonian et al.|2006a I. The decision tree-based machine 
learning algorithm BDT returns a continuous variable called ^ 
which was used to select y-ray-like events. Cutting on this pa- 



for spectral and morphological analysis, respectively ( Ohm et al. 
|2009| ). 

Similar to the H.E.S.S. Standard Analysis, two sets of y- 
ray selection cuts have been defined in Ohm et al. ( 2009] l. For 
the production of sky images the (hard-cuts are used. They re- 
quire a minimum intensity of 160 photo electrons (p.e.) in each 
camera image yielding a superior angular resolution of less than 
0. 1° even at large offsets of 2.5° from the telescope pointing po- 
sition. Additionally, 30% more background events are rejected 
resulting in a 10% higher sensitivity compared to the H.E.S.S. 
Standard Analysis. For the spectral analysis a low energy thresh- 
old is desirable for a broad energy coverage and achieved by 
applying the (^ti-cuts with a 60 p.e. cut on the image intensity. 
For the data set under study, this infers an energy threshold of 
450 GeV for spectral analysis and 700 GeV for morphological 
analysis. 

For two-dimensional sky image generation and morphology 
studies, the template background model (RoweH"2003 |Berge 
et al. 2007|l is applied. For this method the CR background is es- 



' The software which was used to analyse the VHE y-ray data pre- 
sented in this work is the H.E.S.S. analysis package (HAP) in version 
10-06-pl07. 



timated in parameter space rather than in angular space. In this 
analysis, the BDT output parameter ^ has been used to define 
signal and background regions. The normalisation a between 
signal and background is calculated as the fraction of all events 
in the FoV falling into the signal regime, excluding source re- 
gions, divided by the number of all events in the FoV in the back- 
ground regime, again excluding all source regions. The system 
acceptance to measure y-ray like and CR-like events drops off ra- 
dially with the distance to the telescope pointing position. Since 
this acceptance is different for both types of events, a correction 
is applied to a. Sky images obtained with the template back- 
ground model agree with sky images generated with the ring 
background method. The ring background method estimates the 
signal-like hadronic CR contribution at each trial position on the 
sky by integrating events in an annulus centred on that position, 
excluding potential source regions. 

Figure [T] shows the VHE y-ray count map of the region 
around Wd 1 and reveals very extended y-ray emission. The 
complex morphology apparent in the sky image and a potential 
multi-source origin of the emission is investigated and discussed 
in detail in Section[3] Since it is not possible to estimate the back- 
ground from the same FoV due to the fact that observations have 
been carried out within regions of VHE y-ray emission, the On- 
Off background estimation method is utilised to extract spec- 
tral information for the whole emission region, indicated as a 
white circle in Figure [T] Here, the CR background is subtracted 
from the source region [On data) using extragalactic observa- 
tions taken without any VHE y-ray signal in the FoV (Off data). 
To ensure similar observational conditions for On and Off data, 
only On-Off pairs of observations are considered that were taken 
at similar zenith angles and within four months of each other, re- 
sulting in a total live time of 20 hours for the On data set used for 
spectral analysis. The absolute normalisation a between On and 
Off data is calculated using the fraction of total events in both 
observations (again, excluding potential source regions). The re- 
flected background method (Ber ge et al.|2007[ ) is used to derive 
spectral information for smaller regions and the full data set as 
discussed in Section [372l and l33] 

All studies presented in this work were cross-checked by a 
second analysis chain which is based on the H.E.S.S. standard 
event reconstruction scheme (Aharonian et al.|2006a ) using the 
Hillas second moment method (Hillas 1985) and an indepen- 
dent calibration of pixel amplitudes and identification of prob- 
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Fig. 1. H.E.S.S. excess map of the region around Wd 1 corrected 
for the camera acceptance, in units of equivalent on-axis VHE 
y-ray events per arcmin^ and obtained with the template back- 
ground method. The image is smoothed with a 2D Gaussian ker- 
nel with a variance of 0. 12° to reduce the effect of statistical fluc- 
tuations. Significance contours between 4 cr and 8 cr are overlaid 
in black, obtained by integrating events within a radius of 0.22 
degrees at each given position. The green star marks the position 
of Wd 1, the white cross the best fit position of the VHE y-ray 
emission and the black dashed line the Galactic plane. The inlay 
in the lower left corner represents the size of a point-like source 
as it would have been seen by H.E.S.S. for this analysis and the 
same smoothing, normalised to the maximum of HESS J1646- 
458. The dotted white circle has a radius of 1.1° and denotes 
the region which was used for the spectral reconstruction of the 
VHE y-ray emission. Note that the bright region in the lower 
right corner is the source HESS J 1640-465 detected during the 
GPS dAharonian et al.|2006c| l. 

lematic or dead pixels in the lACT cameras. Additionally, the 
Model Analysis ( |de Naurois & Rolland||20(j9] l for the selection 
of y-ray-like events has been utilised to cross-check the spectral 
results. All analyses give compatible results. 

3. VHE Results 

3.1. Position 

Figure [T shows a background-subtracted, camera acceptance- 
corrected image of the VHE y-ray counts per arcmin^ of the 
3°x 3° FoV centred on the best fit position of the y-ray ex- 
cess as obtained with the template background method. The ac- 
ceptance correction has been performed using y-ray like back- 
ground events that pass the y-ray selection cuts. The map is 
smoothed with a Gaussian kernel with a variance of 0.12° to re- 
duce the effect of statistical fluctuations and to highlight signif- 
icant morphological features. Significance contours from Act to 
8cr are overlaid after integrating events within a radius of 0.22° 
at each trial source position. This integration radius is matched to 
the RMS of the Gaussian to resample significant features in the 
sky image and is chosen a priori to match the integration radius 
typically used in the GPS analysis for the search of slightly ex- 
tended sources ( |Aharonian et al.|2006c[ l. Given the extended and 
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Fig. 2. H.E.S.S. radial profile relative to the best-fit position 
of the VHE y-ray emission. The dotted vertical line denotes 
the 95% containment radius used to obtain the spectrum shown 
in Figure |4] Note that the region covering the bright source 
HESS J1640^65 (see Figure [TJ has been excluded for the ra- 
dial profile shown in black by omitting a circle segment with 
220° < 4> < 260° for radii of 1.0° < r < 1.4°. The red graph 
displays the radial profile without excluding HESS J1640-465. 



complex morphology of the VHE y-ray emission the position 
obtained from a two-dimensional Gaussian fit convolved with 
the H.E.S.S. PSF to the raw excess count map obtained for ^hard" 
cuts is used to derive an estimate on the centre of gravity of the 
emission. The two-dimensional Gaussian fit gives a best fit po- 
sition of RA 16H6™50^+27^ and Dec -45°49'12"+7' (J2000). 
Within statistical errors the centre of gravity of the VHE y-ray 
emission is consistent with the nominal Wd 1 cluster position of 
RA 16*'47™00.40^ and Dec -45°51'04.9" (J2000). Based on the 
radial profile shown in Figure |2] the 95% containment radius of 
the VHE y-ray emission relative to the best fit position is deter- 
mined to be 1.1°. This radius is used to extract the energy spec- 
trum presented in Section [33] Note that although the sky image 
gives the impression that the region used for spectral reconstruc- 
tion is contaminated by y rays from HESS J1640^65, this is 
mostly an artifact of the smoothing procedure. The real contri- 
bution is less than 10% in a ring between 1.0 and 1.1 degree 
from the best fit position and only 0.8% in the whole spectral 
extraction region. Within the integration region of 1.1° a total of 
2771 + 139 y-ray excess events at a significance level of 20.9 cr 
pre-trials (20.1 cr post-trials) are found. 



3.2. Morphology 

In order to investigate the multi-source hypothesis two emis- 
sion regions A and B (shown in Figure [3] (left)) are considered. 
The radii of 0.35° and 0.25° of region A and B, respectively, 
are chosen according to the widths of the two substructures. A 
one-dimensional slice in the uncorrected excess image along the 
major axis between the two regions has been produced. The fit 
of two separate sources with Gaussian shape results in a of 
2.0 for 4 degrees of freedom with a probability of 74%. The 
probability that the emission is explained by a single Gaussian 
profile or a constant value is found to be rather low at 0.2% and 
0.1%, respectively. An F-test also supports the multi-source hy- 
pothesis, given that the probabilities that the constant or sin- 



4 



H.E.S.S. Collaboration: VHE y-ray emission from the vicinity of Westerlund 1 



300 




-48°00' 



16''50 



16''40" 



Right Ascension (J2000) 




Fig. 3. Left: H.E.S.S. excess map as shown in Figure [l]but for an enlarged region of 4.5° x4. 5°. The blue-grey boxes denote the 
regions used to generate the one-dimensional slices shown on the right. The white circles denote regions A and B which were 
used for spectral reconstruction (Table |2]i. The weak y-ray emission seen in the lower right corner next to HESS J1640^65 is 
HESS J1634-472, also detected during the GPS ( Aharonian et al.|2006c[ l. Right, top: Distribution of VHE y-ray excess events in 
the blue-grey box, oriented along the two emission regions associated to HESS J1646^58 and starting at low Declination angles. 
The results of a fit of a constant and of two sources with Gaussian shape are indicated as black and red line, respectively. Right, 
bottom: Same as top, but oriented along the Galactic plane, starting at low Right Ascension angles, close to the bright source 
HESS J1640^65 at li = -1.25°. The slice has been truncated at 350 excess events in order to emphasize the VHE y-ray emission 
from HESS J1646-458. 



gle Gaussian emission models are preferred over the double 
Gaussian fit are < 0.02 and < 0.01, respectively. 

Figure[3](left) also suggests a contribution from diffuse VHE 
y-ray emission along the Galactic plane which extends 1° to 
2° from region A north-eastwards. This impression is supported 
by the one-dimensional slice shown in Figure [3] (bottom right), 
where the significance of the emission in all bins with distance 
0.5° < d < 1.8° from the centre of region A is between 2 cr and 
Act. This difiiise emission could be due to unresolved VHE y- 
ray sources or a Galactic diffuse emission component, caused by 
the interaction of GCRs with molecular material located along 
the Galactic plane. As will be shown later in Figure 6] there is in- 
deed molecular material located in this region which could act as 
target for the interaction with CRs (as described in e.g. 'Casanova' 
et al. 2009 |l and could account for part of the observed emission. 



The statistics of region A and B compared to the entire source re- 
gion as obtained with the template background method are given 
in Table [U 

The studies presented here show some evidence for a multi- 
source morphology and a separation into multiple VHE y-ray 
sources. Moreover, spectral variations across the whole emis- 



sion region, e.g. as observed for HESS J1825-137 (Aharonian 
|eraL]|2006d| ) or in the case of the H.E.S.S. sources in the 
Westerlund 2 field ( [Abramowski et aL||2011| l, could be appar- 
ent, which would further support the multi-source hypothesis. In 
this case, an energy-dependent morphology can be expected. To 
test this hypothesis, the complete data set has been divided into a 
low-energy band, containing events with reconstructed energies 
£ < 1 .0 TeV and a high-energy band, containing events with re- 
constructed energies £ > 1 .0 TeV. The unsmoothed excess maps 



in coarse bins of 0.3° width are used to test the underlying distri- 
bution. Ax^ test is performed using the number of excess events 
in each bin and reveals a value of 95.4 for 76 degrees of free- 
dom (four bins covering HESS J1640^65 have been excluded 
in the calculation). Prior to the test a. probability po of 0.05 
is defined to accept/reject the null hypothesis. The p-value of 
the test is 6.5% > po, such that the null hypothesis that both 
excess maps follow the same underlying distribution cannot be 
rejected. Although no energy-dependent morphology can be de- 
tected from this test, a multi-source origin is preferred given the 
low probabilities for the single source fits. 

3.3. Spectrum 

The spectrum obtained for the whole emission region us- 
ing the On-Off background estimation method is shown in 
Figure |4] In the fit range between 0.45 TeV and 75 TeV, the 
spectrum is well described by a power law: dN/dE - Oo x 
(E/l TeV)-r with a photon index F = 2.19 ± 0.08,tat + 0.20,sys, 
and a differential flux normalisation at 1 TeV of Oq = (9.0 + 
1.4stat + 1.8sys) X lO^'^TeV^'cm^^g-i j^^^^. translates into 
an integral flux above 0.2 TeV of F(> 0.2 TeV) = (5.2 + 
0.9) X 10"" cm"2 s'K The x^ for the power law fit is 9.9 for 
7 degrees of freedom, yielding ax^ probability of 19%. 

Additionally, the differential energy spectra for region A 
and B are determined using the reflected background estimation 
method with results found to be consistent with the On-Ojf back- 
ground technique. The integral flux above 0.2 TeV as well as the 
spectral results obtained from a power-law fit are summarised 
in Table |2] and compared to the results for the spectral analy- 
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Table 1. VHE y-ray statistics for the regions shown in Figure fllandls] 



Region 


RA (J2000) 


Dec (J2000) 


e 


On 


Off 


a 


Excess 


Significance 




deg 


deg 


deg 


events 


events 




events 


a 


Full 


251.856 


-45.909 


1.1 


19032 


1107471 


0.014682 


2771 


20.1 


A 


251.370 


-45.585 


0.35 


2313 


120104 


0.014998 


511 


10.0 


B 


251.682 


-46.513 


0.25 


1149 


58995 


0.014876 


271 


6.5 




Fig. 4. Top: Differential VHE y-ray energy spectrum of 
HESS J1646^58. The data are fit by a power law dNIdE = 
<E)o X (£■/! TeV)""^. Arrows indicate the 95% upper limits for 
spectral bins which are compatible with a zero flux within Icr. 
Also shown are spectra as obtained for region A and B. Bottom: 
Residuals of the power-law fit for the 1.1° y-ray emission region. 



sis of the whole 1.1° region. The differential energy spectra for 
these two regions are shown in Figure |4] as well. Within statis- 
tical errors, no change in photon index between the three stud- 
ied regions is apparent, further supporting the lack of energy- 
dependent morphology across the source based on the current 
data. 



4. Discussion 

In this section, the spectral and morphological results are 
used to elaborate possible acceleration scenarios related to 
HESS J1646^58. Although the morphological analysis prefers 
a two-source approach (at the ~ 2.5 cr level), the similarity of 
the spectra from region A &lB does not allow to substantiate the 
preference for a two-source scenario further Owing to this am- 
biguity, we pursue both possibilities through all portrayed coun- 
terpart scenarios. 

Due to the large size of HESS J1646-458 (~ 2°), a few ob- 
jects that could be regarded as potential VHE y-ray emitter are 
found in the region of interest. Besides the massive stellar clus- 
ter Westerlund 1, a high spin-down power pulsar PSR J 1648- 
4611 ( Manchester et al.||2005[ l which has recently been discov- 
ered to be a high-energy y-ray pulsar with a possible unpulsed 
y-ray component is found ( Kerr 2009] Abdo et al.|2010 2011 



Fig. 5. H.E.S.S. smoothed VHE y-ray excess contours in blue 
overlaid on the Molonglo 843 MHz map (grey scale, in Jy/beam; 
Green et al.|1999| l. Also shown are S NRs (purple, , Green 2009) , 
the pulsar PS R J1648-461 1 (M anches ter et al.|2005| l, the LMXB 
4U 1642-45 dForman et aL||1978 ), Fermi-LAT sources (Abdo 



etaD[20TT) the magnetar CXOU J164710.2-455216 (Muno 
et al.||20()6a i as a black upright triangle and the stellar cluster 
Wdl. 



455216 ( |Muno et al.| |200"6al l and three unidentified Fermi- 
LAT sources, 2FGL J1 650.6-4603c, 2FGL J1651.8-4439c and 
2FGL J1653.9^627c jAbdo et al.|r2011 ) are located within 
HESS J 1646-458 (Fig.|5f~^ 

In the following, the discussion will focus only on the astro- 
physically associated objects with known distance and energet- 



ics. 



4.1. 4U 1642-45 



The X-ray bright LMXB 4U 1642-45 with an accreting neu- 
tron star lies within the strong emission r egion A (Fig. [5]l and is 
located at a distance of 8. 5 to 11.8 kpc (|van Paradijs & White 
[T9951 IGrimm et al.|[2002t [Gilfanov et al.||2003| l. It exhibits an 
average X-ray luminosity ofLx ~ lO-^'^ergs"^ and Grimm et al. 
( 2002[ l reported of periods of super- Eddington luminosity from 
ASM observations. 

So far, only high-mass X-ray binary systems have been as- 
sociated with point-like and variable VHE y-ray sources (e.g. 
LS 5039, [Aharonian et al.|[2005| |2006e[). At 8.5 to 11.8 kpc. 



Furthermore, the low -mass X-ray binary (LMXB) 4U 1642-45 
( |Forman et al.|1978| l, the magnetar candidate CXOU J 1647 10.2- 



HESS J1646^58 would be 320 to 450 pc in size and region 
A 100 to 144 pc, for angular diameters of 2.2° and 0.7° respec- 
tively. These inferred source sizes are well beyond the H.E.S.S. 
PSF in this analysis (26 to 35 pc in diameter). A check for vari- 
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Table 2. Spectral properties of the different TeV extraction regions. 



Region 




Oo(lTeV) 




r 


F(> 0.2 TeV) 


% total 




10- 


'2 TeV-' cm-2 s"' 






10-"cm-2 S-' 




Full 


9.0 


± 1.4 


2.19 


±0.08 


5.2 ±0.9 


100 


A 


2.1 


±0.3 


2.11 


±0.12 


1.1 ±0.2 


21 ±4 


B 


1.4 


±0.2 


2.29 


±0.17 


0.8 ± 0.2 


15 ±4 



ability at the position of 4U 1642^5 did not reveal any indica- 
tion of such as has been observed, e.g. in the case of LS 5039 
and HESS J0632+057 (Aharonia n et al.|2006e||2007 ^. 

In summary, an association of 4U 1642-45 and the entire 
VHE y-ray emission region is unlikely. The extent of the subre- 
gion A also disfavours a scenario with 4U 1642-45 and a sub- 
region as well. An association with region A would be solely 
based on spatial coincidence but would then present a new class 
of VHE y-ray sources. 

4.2. CXOU J1 6471 0.2-45521 6 

This anomalous X-ray pulsar is considered a magnetar and 
has been associated with Wd 1 given its apparent proximity 
to the cluster and the low probability of a random association 
( Muno et al.||2006^ . Shortly after its X-ray discovery, an out- 
burst in X-rays was reported from Swift observations (Campana] 
& Israel|[2006| l. Additionally, XMM-Newton observations were 



conducted prior to and after the outburst (Muno et al. 2007 
[Israel et al.|2007l ). 

The observed X-ray luminosity is Lx ~ 3x10^^ erg s"' which 
increased during the outburst by a factor of ~ 100. The rota- 
tion period is 10.6s with a derivative of P = 9.2 x 10 '^ss ' 
which infers an X-ray spin-down power of £ = 3 x 10^' erg s ' 
and a characteristic age of t = 1.8 x lO^yrs. The surface 
magnetic-field strength is ~ 1 x lO'^^G. The relevant properties 
of CXOU J 1647 10.2-4552 16 are listed in Tabled 

Since the X-ray luminosity exceeds the rotational spin-down 
power, this suggests that the observed X-ray emission is not 
due to spin-down processes, but some other form. For magne- 
tars, the intense magnetic field can yield considerable power 
in magnetic dissipative processes. Following Zhang ( 2003[ l, 
such power can be of the order L^ag ^ 10^^ erg s~', which, 
for CXOU J164710.2^55216, would" dominate over the rota- 
tional spin-down power If acceleration of particles to TeV en- 
ergies occurs due to magnetic energy release one could expect a 
rather compact TeV-emission region similar to spin-down pow- 
ered TeV pulsar wind nebulae (PWNe) like the Crab nebula 
( [Aharoniari et al.|2006a| l. 

However, it is still under investigation if magnetars can 
exhibit PWNe in general. For most pulsars, the respective 



PWNe would be too faint to be detected in X-rays ( Gaensler 



& Slane 2006 1. In the VHE y-ray regime, young (0(10^ yrs)) 
and high spin-down power {0(10^'^ erg s ')) magnetars seem to 



be promising objects ( Halpern & Gotthelf|2010 1 

However, the rotational spin-down power of 
CXOU J164710.2-455216 is too low to account for both 
the entire observed VHE y-ray emission and either subre- 
gion. Therefore, any acceleration of particles would have to 
involve a magnetic energy release with a power output of 



-^mag 



10^5 erg s-1. 



The observed emission could then be 



explained by an energy conversion process that operates with 
efficiency up to 100% but the size of HESS J1646-458 (-160 pc 
at the distance of the magnetar) stands in contradiction to the 
expected compact region. Either of the two subregions could be 



energetically explained with efficiencies in the order of 0.5 but 
again as one would expect the resulting PWN to be compact and 
close to the magnetar itself, these scenarios are disfavoured as 
well. 

All in all, VHE y rays from CXOU J164710.2^55216 is 
not a favoured scenario to account for HESS J1646-458 or parts 
from it as its relevant properties and current (V)HE y-ray obser- 
vations do not support such an approach. 

4.3. PSRJ1648-4611 

Among known VHE y-ray sources, PWNe are the most abundant 
source type: Roughly one third of these sources are associated 
with PWN systems (e.g. Hinton & Hofmann 2009| l. The pulsar is 
found to be located in the centre of the nebula (e.g. the Crab neb- 
ula, ^^^omanetaiT'2006a') or offset from it (e.g. the Vela X neb- 
ula, Aharonian et al. 2006b). Most of the VHE y-ray-emitting 



PWNe are spatially extended and offset from the pulsar position 
and efficient in terms of converting available spin-down power 
into y-ray emission. In general, PWNe firmly associated with 
known pulsars convert 10% at most of their spin-down power 
into y-ray luminosity ( Gallant^2007^) . 

The high spin-down power pulsar PSR J 1648- 46 11 with 



10^^ ergs ' ( [Manchester et al" 2005 i is located 



£■ = 2.1 X 

within the emission region HESS J1646^58 and associated with 
the Fermi-LAT source 2FGL J1648.4^612 ( |Abdo e t al. 2011). 
Pulsed y-ray emission with the rotation period of PSR J 1648- 
4611 as well as signatures for a constant (possible PWN) con- 
tribution have been reported. Therefore, this object has to be 
considered a y-ray pulsar at low GeV energies. Inspection 
of archival Chandra data with lOks exposure (obsid: 11836) 
shows no indication for the presence of an X-ray counterpart 
to PSR J1648^611. Following the relation between spin-down 
power and X-ray luminosity of PWNe ( Mattana et al.|2009 ), the 
expected X-ray luminosity Lx ~ 5x10-^" erg s~' of the associated 
PWN would not be detectable with current X-ray instruments 
and is hence consistent with the non-detection with Chandra. 

The first attempt to explain HESS J1646-458 is a single- 
source scenario which is motivated by the presence of the HE 
pulsar PSR J1648-4611. The inferred y-ray luminosity (0.1 - 
100 TeV) at the distance of PSR 11648^611 would require a 
conversion efficiency eeff ~ 1. Additionally, the size of the 
VHE y-ray-emitting region would extend over roughly 200 pc 
which is a factor 3-10 larger than known extended VHE 
y-ray-emitting PWN systems. Given the size, cooling losses 
would lead to considerable softening of the VHE y-ray spec- 
trum and a downward shift of the maximum energy with in- 
creasing distance from the pulsar - similar to the spectral soft- 
ening of HESS J1825-137 (Aha ronian et al. 2006f). The VHE 
y-ray photon index reconstructed at the position of PSR J1648- 
4611 is F = 2.37 + 0.43 using the re^ecfetZ-background regions 
method and assuming a point-like source origin. This is com- 
patible within statistical errors with the emission of the entire 
region. However, the available statistics do not permit firm con- 
clusions about potential spectral changes across HESS J1646- 
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Table 3. Properties and inferred y-ray luminosities of the pulsars and the VHE y-ray emission regions. 



Object / Region 


T 


d 


Po 


P 


E 






10' yrs 


kpc 


s 


lO-'^ss-i 


lO^-* ergs-' 


lO^-* ergs-' 


CXOU J164710.2^55216'' 


1.8 


5 


10.6 


9.2 


3 X 10-' 


26 


PSR J 1648-46 11* 


1.1 


5.7 


0.2 


0.2 


21 


34 


HESS J1646-458 




4.3 








19 


Region A 




4.3 








4.6 


Region B 




4.3 








2.9 



Notes. Obtained by scaling the observed VHE y-ray flux between 0. 1 and 100 TeV to the nominal or adopted distance of the respective object. 
References. '°' |Muno et al.|j2006a| >; [lsrael et al.|j2007) . [Manchester et al.| ( |2005 i. 



458 to be drawn. In summary, the unprecedented high efficiency 
needed and the size of the VHE y-ray emission region disfavour 
HESS J1646-458 as a very extended PWN. 

In a two-source approach as motivated in Section |3.2| a 
displacement of PSR 11648^611 from either subregion is ap- 
parent. The pulsar is displaced by ~ 50 - 70 pc from region 
A and B. These offsets are large compared to those of known 
VHE y-ray PWNe but could in principle be explained by rela- 
tive proper motion of the pulsar assuming a transversal veloc- 
ity of C)(500kms-') which is at the upper end of the range of 
known transversal motions of pulsars ( [Manchester et al.|2005[ l. 
In this scenario, one of the two subregions could in principle 
be powered by the pulsar with egff a; 0.1 . . . 0.2. In this case, the 
morphology would reflect the ambient conditions, e.g. one of the 
subregions could be the result of the expansion of the PWN into 
a low-density medium or could be due to an asymmetric reverse 
shock of the SNR. 

Recently, |Luna et al.| ([2010) proposed that a region of low 
density in CO seems to partially match structures in prelimi- 
nary VHE y-ray data (Ohm et aT]|2010b| l and for which a SN 
event with PSR J 1648-46 11 as precursor could be responsible. 
However, the kinematic age of the cavity is about 55 times larger 
than the characteristic age of PSR J 1648-461 1 and the inferred 
subsonic expansion velocity is insufficient to accelerate particles 
up to the VHE y-ray regime. Moreover, other SNR candidates at 
the position of this cavity are not to be found in archival data and 
the morphology of HESS 11646^58 does not strongly motivate 
a shell-like structure centred at about 16M7"'23 .3% -45°42'5.2" 
with an inner radius of ~0.5° (|Luna etal.|2010 1. Although it can- 



not be ruled out that this cavity could be blown as a consequence 
of a SN it seems unplausible that the whole ~ 25 pc-thick shell 
would expand uniformly at 6 - 8kms-' and thereby giving rise 
to particle acceleration up to TeV energies. 

Similarly as concluded for the magnetar, HESS J 1646^58 
seems unlikely to be explained as a very extended PWN pow- 
ered by PSR J 1648-461 1. Either one of the subregions could be 
explained by an offset PWN. The infeiTed offsets would be com- 
paratively large but the required efficiency would be amongst 
known TeV PWNe. 



4.4. Westerlund 1 

The motivation for H.E.S.S. observations of the Wd 1 region 
has been outlined earlier (see Section [T} and the analysis of 
HESS J1646^58 revealed that Wd 1 is located close to the cen- 
troid of the VHE y-ray emission. Hence this massive cluster is 
an interesting object to consider. As some of its properties are 
still subject of on-going investigation certain assumptions are 
made: (1) Due to the high visual extinction of A,, » 12™^^^ 
distance estimates based on photometry and spectroscopy vary 



o 
o 
o 



o 

u 
Q 



-45°00'- 



-46°00'- 



-47°00 




16''55' 



16"50' 



16"45'" 16"40"' 
Right Ascension (J2000) 



Fig. 6. HI 21cm line emission map at a radia l velocity of 
about -55kms-' ( [McClure-Grifliths et"aL]|2005| between 20 
and 80 K. Bright regions are Hl-depleted whereas dark re- 
gions are Hl-dense. Overlaid are the CO contours (red) at 12.5, 
22.5 and 32.5 K (Dame et al. |2001 1 along with the smoothed 
HESS 11646^58 contours in blue. The estimated position of the 
HI void B3 is marked by the small green dashed circle (Kothesj 
|& Dougherty|2007] ). The large green dashed circle indicates the 
region used for the spectral reconstruction in the TeV regime and 
to compute the ambient matter density. The position of Wd 1 is 
marked by a white star in the centre. The circles in the lower left 
comer indicate the beam sizes of the respective radio observa- 
tions. 



strongly between 1.1 kpc and 5.5 kpc ( Westerlund'l l987[ Piatti 



et ai:] [T998l [Clark et al . 2005; Br andner e t al. 2008 ). IKoflieFfe 



[Dougherty (2007 ) derived a distance of 3.9 kpc based on the Hi 
line emission which , Luna et al.|(2010) extrapolated to 4.3 kpc 
using the lAU distance recommendation to the Galactic Centre 
of 8.5 kpc. In this work, we adopt the 4.3 kpc albeit using a more 
recent Galactocentric distance estimate of 8.33 kpc ( [Gillessen 



|et al. [2()09| . (2) Different model-dependent approaches to es- 



timate the age suffer from the apparent lack of main-sequence 
stars. Recently, Negueruela et al.[ (2010l estimated the age to be 
> 5Myrs and an age of 5Myrs is adopted here. (3) A cluster 
mass of 10'' is assumed. 

Estimates based on the stellar population of Wd 1 imply that 
80 to 150 stars with initial masses exceeding 50 Mq have already 
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undergone a SN explosion ( |Muno et al.||2006bl l. However, nei- 
ther radio nor X-ray observations have revealed any SNR candi- 
date. In order to estimate the total kinetic energy dissipated by 
the system through SNe and stellar winds the Starburst99 clus- 
ter evolution model ( Leitherer et al.|1999 Vazquez & Leitherer 



2005 Leitherer et al.|2010| l has been used. The default parame- 
ters with the standard Kroupa initial mass function (with the ex- 
ponents 1.3 and 2.3, (Kroupa 2002 )) and the evolutionary-model 
option, have been chosen. The total energy dissipated at the nom- 
inal age of Wd 1 including stellar winds and SNe amounts to 
£kin = 3.0x10^3 erg. 



For an ad iabatically expanding wind (Weaver et al. 1977 
SiUch et al. 2005) , a structured and expanding (~ 30km s~') hot 
bubble with an outer shock radius of ~ 250 pc is expected to 
form around Wd 1 even if only a moderate number of 50 OB 
stars are considered. However, a dedicated search for such a su- 
per bubble-like structure in HI and CO data at the radial ve- 
locity of Wd 1 did not reveal any signatures of such a feature 



( [McClure-Griffiths et al.|2005||Kothes & DoughertylSOOT] !. The 
latter authors find indications for a much sm aller (55 pc) and 
slowly expanding (3kms ') feature (B3 in Kothes & Dougherty 
2007| l, see Fig. |6] However, in this first approximation, radiative 
cooling is not considered although this cooling process might 
affect the evolution of stellar cluster winds (e.g. discussed in 

IWiinsch et airSOOSjl. 

[Kothes & Dougherty ( 2Q07\ interpret the formation of the 
HII region complex G340.2-0.2 (depicted in Fig. 6} as triggered 
by the stellar wind activity of Wd 1 . Indeed, if Wd 1 is located at 
a distance of 4.3 kpc some correlation between VHE y rays and 
the location of this HII region is expected. 

An average gas density for atomic hydrogen nei and molec- 
ular gas can be derived using available HI and CO data in 
the velocity range of Wd 1 (-50 to -60 km s"') for the entire 
1.1° region of interest. For this, all pixel values in this velocity 
range and within 1.1° from the cluster position are considered as 
well as the SGPS beam size (130arcsecs) and the oversampling 
factor (11.97). With the HI int ensity-mass conversion factor of 
1.823 X 10''^cm-V(Kkms-i) (Yamamoto et al. "2003), the re- 
sulj^is hhi = 0.24 cm"3. In the CO data ( Dame et al. 2001) , a 
low-density region around the cluster seems to be apparent as 
well (Fig. [6]l. This feature is larger than B3. With a beam size 
of 450 arcsecs and a oversampling factor of 1, the application of 



the conversion factor of CO to H2 of 1 .5 x 10^° cm ^/(Kkm s ') 



,-2 / 



( Strong et al.i2004 ) leads to an average density of H2 molecules 
in units of atomic hydrogen of hh, = 12.16cm 3. Note that the 
X factor used for the CO data incorporates caveats pointed out 
Strong et al. ( 2004[ l. The required energy in CRs to power the 



m 



y-ray emission can now be estimated: 



£cR = 2.1 X 10^' 



IGeV 



5.8 X 10^5 erg s- 



(124^) 



where L>i cev is the high-energy luminosity computed with the 
spectral results presented earlier between 1 GeV and 1 PeV. 

When interpreted in a single-source scenario, the required ef- 
ficiency for transferring kinetic energy in shocks or turbulences 
into energetic particles through acceleration is therefore at the 
level of ~ 0.1%. For a similar argument as for the PWN interpre- 
tation, a leptonic origin of the VHE y-ray emission is difficult to 



^ The density n is proportional to 2 b'PiR^ fs^ where p, are the pixel 
values of the HI or CO map, R the radius of the region of interest, b the 
beam size of the respective experiment and is the oversampling factor 
in the respective HI or CO maps used. 



reconcile with the large extent of HESS J1646^58 which trans- 
lates into a size of 160pc at the distance of the cluster Given the 
large photon density in the environment of Wd 1, a fast convec- 
tive transport of the electrons would be required to prevent them 
from cooling. Hence a dominant hadronic origin is favoured in 
this approach. 

A possible two-source scenario would also consist of a dom- 
inant hadronic CR component as the stellar photon field would 
lead to a rapid cooling of VHE electrons. In this case, bright 
VHE y-ray structures would trace dense features in HI and CO. 
In particular, region A lies at the edge of B3 and is coincident 
with dense structures in HI (see Fig. |6]l which would naturally 
provide sufficient target material. Region B, however, remains 
comparatively dark as HI and CO data do not suggest a higher 
abundance of target material. 

In summary, Wd 1 and its massive stars favour a hadronic 
mechanism for the entire emission region. Here, the size and the 
inferred energetics of HESS J1646-458 could be plausibly ex- 
plained. In case of the two subregions, or multiple source regions 
in general, more observational data in all wavelength bands and 
detailed modelling are required. 



5. Summary and Conclusions 

In summary, HESS J 1646-45 8 is a new VHE y-ray source 
found towards the unique massive stellar cluster Westerlund 1 
and a number of other potential counterparts. The large size of 
HESS 11646^58 however, over 2 degrees in diameter making 
it one of the largest TeV sources so far detected by H.E.S.S., 
presents a challenge in identifying a clear counterpart (or a num- 
ber of counterparts) to explain the VHE y-ray emission. 

The detection of degree-scale VHE y-ray emission, namely 
HESS J1646^58, towards the stellar cluster Westerlund 1 with 
a total significance of 20 cr from H.E.S.S. observations per- 
formed in the years 2004, 2007 and 2008 (33.8 hrs live time) 
is reported. The energy spectrum between 0.45 TeV and 75 TeV 
of the entire region is best fit by a simple power law with 
an index F - 2.19 + 0.08stat ± 0.2sys and a normalisation at 
ITeV Oo = (9.0 ± 1.4,s,at + 1.8sys) x lO^'^ TeV-'cm-^g-' with 
;^'^/ndf = 9.9/7. The integrated flux above 0.2TeV amounts to 
(3.49 + 0.52) X 10 " cm ^s-i. The VHE y-ray luminosity be- 
tween 0.1 and 100 TeV is 1.9 x 10^^ (t//4.3kpc)2ergs-i. 

The centroid of HESS 11646^58 is consistent with the nom- 
inal position of Wd 1. The observed VHE y-ray emission region 
has a diameter of about 2° which translates into a spatial ex- 
tent of 160 pc at the distance of Wd 1 or 200 pc at the distance 
of PSR J1648-4611. hi either case, the size of HESS J1646- 
458 would be the largest among currently known VHE y-ray 
sources, if the emission is of a single-source origin. This is sup- 
ported by the lack of spectral changes across HESS J1646-458 
within statistical errors. Although there is some evidence for a 
multi-source morphology, the limited statistics hamper a detailed 
investigation into the presence of multiple sources. 

In a scenario where one astrophysical object is responsible 
for HESS 11646^58, Wd 1 could naturally account for the re- 
quired injection power provided that about 0.1% («/12.4cm"-') 
of the kinetic energy released by stellar winds and supernova ex- 
plosions are converted into particle acceleration. In this case, the 
stellar wind and SNe activity of Wd 1 would strongly affect the 
surroundings but also reflect the ambient conditions. 

In a split of HESS 11646^58 into two distinct subregions, 
however, a superposition of two or more sources adding up to 
the observed VHE y-ray emission region could be possible. In 
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this case, one of the subregions could also be explained by a 
PWN with comparatively large offset from PSR 11648^611. 
Despite the spatial coincidence, the LMXB 4U 1642^5 and the 
magnetar CXOU 164710.2-455216 are not likely to account for 
HESS J1646-458 or parts of it. 

Further multiwavelength observations (in particular those 
covering the VHE y-ray source in full), and deeper VHE y-ray 
coverage will no doubt be valuable in shedding light on the na- 
ture of this source. 
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